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Extending the scope of microscopic solvability: Combination of 
the Kruskal-Segur method with Zauderer decomposition 
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PACS 47.54.-r - Pattern selection; pattern formation 

PACS 11 . 10 . Jj - Asymptotic problems and properties 

PACS 81 . lO.Aj - Theory and models of crystal growth 

PACS 68.70.+W - Whiskers and dendrites (growth, structure, and nonelectronic properties) 

Abstract. - Successful applications of the Kruskal-Segur approach to interfacial pattern forma- 
tion have remained limited due to the necessity of an integral formulation of the problem. This 
excludes nonlinear bulk equations, rendering convection intractable. Combining the method with 
Zauderer's asymptotic decomposition scheme, we are able to strongly extend its scope of appli- 
cability and solve selection problems based on free boundary formulations in terms of partial 
differential equations alone. To demonstrate the technique, we give the first analytic solution of 
the problem of velocity selection for dendritic growth in a forced potential flow. 



The fundamental equations describing the growth of a 
crystal into its undercooled melt are very difficult to solve, 
if surface tension effects are accounted for, even when re- 
striced to the simplest case of merely diffusive heat trans- 
port. On the other hand, the capillary length do describ- 
ing these effects is typically very small in comparison with 
other length scales of the problem such as the sizes of grow- 
ing patterns or the diffusion length. Therefore, it was a 
natural step to first look for solutions with do set equal 
to zero. This simplified problem was solved exactly by 
Ivantsov [1] who showed that the crystal can grow in the 
shape of a parabola in 2D or a paraboloid in 3D. A ma- 
jor drawback of these solutions is that they constitute a 
whole continuum for any given undercooling: the mathe- 
matics fixes only the Peclet number P c — Vp/D, where V 
is the growth velocity of the crystal, p the tip radius of the 
parabolic needle, and D the thermal diffusion coefficient. 
Hence, only the product of velocity and length scale is de- 
termined, but neither of the two quantities separately. In 
experiments, a given undercooling leads to both a well- 
defined growth velocity and a well-defined tip radius of 
the needle crystal, which after developing side branches is 
called a dendrite. This situation became known as the se- 
lection problem of diffusion-limited dendritic growth and 
is was not solved until some twenty years ago [2-5], with 
the advent of microscopic solvability theory. 

Because the theory was mathematically complex and 



not very intuitive, it failed to enjoy unanimous appraisal. 
Moreover, its success in explaining experiments remained 
controversial to some extent [6]. It has been emphasized 
by Tanveer [7] that even small fluid flows in the melt might 
account for changes in the theoretically predicted scalings 
as the problem is structurally unstable. Hence, selection 
theory should be extended to nondiffusive transport such 
as convection. To our knowledge, the only approach to 
solvability theory available so far for models with convec- 
tion is due to Bouissou and Pelce (BP) [8]. Their method 
relies on a linearized solvability condition, which prevents 
it from becoming exact in the limit of vanishing do. Also, 
it has been shown [9] that nonlinearity may be crucial in 
problems involving multiple parameters. Hence, a method 
would be more than desirable that takes nonlinear solv- 
ability into account but can also deal with problems not 
permitting an integral formulation. This letter presents 
such an approach. 

Lengths are nondimensionalized by the tip radius p of 
the Ivantsov parabola [1], which is then given by y — yo = 
i — ix 2 in a comoving coordinate system. We use con- 
formal parabolic coordinates x = £17, y — yo — \{i] 2 — £ 2 ), 
so the Ivantsov parabola reads 77 = 1. The relevant 
nondimensional parameters of the problem are the growth 
Peclet number P c , defined above, the stability parame- 
ter cr, and the flow Peclet number Pf. The latter are 
given in terms of physical quantities by a = 2doD/[Vp 2 ], 
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Pf = Up/D, where U is the velocity of an imposed flow. 

In order to demonstrate the power of the method, we 
first show how it simplifies a solved problem, the growth of 
of a needle crystal under purely diffusive heat transport. 
Next we deal with the hitherto unsolved selection problem 
of a crystal growing in a potential flow, where the basic 
field equations are nonlinear. We simplify the presenta- 
tion using some approximations that can in principle be 
avoided and will be removed in a more extended version of 
this article. The first problem was treated via the Kruskal- 
Segur method [10] by Ben Amar et al. [3]. They start from 
an integral equation describing steady state growth in the 
limit P c — > 0. The integral term is linearized about the 
Ivantsov parabola, setting r] s = 1 + h [%(£) is the inter- 
face position]. After reducing the integral to a local ex- 
pression using sophisticated complex analysis [5] , one finds 
the dominant behavior of the solution near a singularity 
in the complex plane at £ = — i: 



uk = (i - 
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where 



k = 



(e + (i + m 2 ) 
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(l + /l' 2 )5 
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is the curvature [the prime denotes a derivative w.r.t. the 
argument (i.e. £)]. Equation (1), a second-order nonlin- 
ear differential equation for /i(£), contains all the infor- 
mation needed to compute the transcendental corrections 
(i.e., the mismatch function) that have to be suppressed 
at the needle tip for selection to be possible. Essential for 
its derivation was the use of an integral equation, available 
only for linear bulk equations. 

We now rederive Eq. (1) from the differential equa- 
tion formulation of the free-boundary problem directly. 
The field equation for the problem linearized about the 
Ivantsov solution is just the Laplace equation d^T+d^T = 
0, both in the liquid and solid phases. The interface 
boundary conditions become: 
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T\i = T\ s + h, 



h + ih 1 + (d 71 -h'd^){T\ l -T\ s ) = o. 



(3) 



Of these, the first equation is the Gibbs-Thomson relation, 
the second expresses continuity of the temperature at the 
interface (subscripts I and s refer to evaluation at the po- 
sition of the Ivantsov parabola) , the third is the continuity 
equation. Noting that (d| + <9 2 ) = (d$ +id v )(d^ —id n ), we 
replace the bulk equation with 



(d^ + id v )T = in the liquid 
(d i - id v )T = in the solid . 



(4) 
(5) 



From Eqs. (3), (4) and (5) we get, neglecting terms 
quadratic in h! [11] 



and after one integration, obtaining the integration con- 
stant from the boundary condition h — ► (£ — ► oo), we 
recover (1) almost effortless. 

Let us briefly discuss the philosophy of this approach. 
The temperature field satisfies the Laplace equation, 
solved by T(£, rf) = h [£ + Hv ~ 1)] + h K - Hv ~ !)]■ After 
inserting this general solution into the boundary condi- 
tions (3), valid at r\ = 1, we analytically continue these to 
the vicinity of w = — i (w is the analytic continuation of £). 
Some of the terms must become singular there to compen- 
sate for the singularity of the curvature term. Since the 
solution for the liquid must be analytic in the upper half 
plane (77 > 1), the f 2 term remains regular near w = —i 
and is hence negligible. The other term is a solution of 
(4). Similarly, dropping regular terms from the solid-side 
solution, we keep the term that solves (5). This procedure 
gives a valid approximation near the singularity. Far away 
from singularities it is also justified, because the curvature 
term in Eqs. (3) can be linearized and the corresponding 
inhomogeneity ignored. We also have to take care of the 
singularity at w = i. Since the final result has to be real, 
this singularity leads to the complex conjugate. 

Let us determine the transcendental mismatch for later 
reference when considering the case with flow. Far from 
the singularity, the homogeneous part of the linearized 
Eq. (1) reads 

a ( h " + TTe h ') + v/TT ^ (1 ~ ® h = • (7) 

This may be solved using WKB techniques, yielding an 
outer WKB solution: 

h = B(l-iO~*(l + iO _i exp{a-iS(0},(8) 
S(0 = i J{l-i?)i(l + i?)id$'. (9) 

—i 

To obtain the inner equation near £ = —i, we rescale 

h = a a 4>, £=-i + i<7 a t. (10) 
From equation (1), we have 

a 2a t<f> = ^-^K + O (a-3 a+1 ) , (11) 

where 

~ (2t + 20)3 { (1 - 2 )i + (2i + 20)(1 - 2 )2 J ' 

(12) 

To balance both sides of equation (11), we need a = |, 
hence t<f> = K, and the asymptotic behavior of <j> for large t 
is 4> ~ 2~2t~f . Linearizing about the asymptotic solution 
and performing a local asymptotic analysis near t — 00, we 
get the decreasing eigenfunction g = t - i cxp j — |2*i* j, 
leading to the transcendental behavior 



(1 -«£)/» = -2 (T\s)' = an' , 



(6) 



= 2~H~i +Ag(t) (t->oo) 



(13) 
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where the nonlinear eigenvalue A can be found by solving where w = (Tj, T v , tp) and 
t<j) = K numerically [3], imposing the calculated asymp- 
totic behaviour. The constant B may be related to the 
nonlinear eigenvalue A by matching the outer WKB cor- 
rection (8) with the transcendental contribution found in 
the inner domain (13): A = 2~i a~^s B. For the tip slope, 
we get 
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Using the Kronecker symbol Sij, the remaining matrices 
may be written Bij = 6n [~Pf£,dj\ + Pf(rj — l)5j2 



= 2T (T -Mi m (A)exp{ ( T-5 5'(o)| , (14) -e ^ Pf(ri 1)2 5nS j3 , on, = 6^(1 - 6 a ), # 



dj — 



a result that shows that with isotropic capillary length 
there is no solution to the selection problem, as the right- 
hand side of (14) is different from zero. 

Let us now consider how the approach works when a 
potential flow is imposed externally, with the nondimen- 
sional flow velocity U tending to ~Pfe y for y — > oo. We 
choose this irrotational frictionless flow, because for this 
case the exact solution of the problem without surface ten- 
sion, i.e., the analog of the Ivantsov solution, is known [12]. 
Clearly this is just a toy model, as it cannot be used to 
study viscosity effects. But it is a useful simple example 
in demonstrating our method, extensible to more realistic 
flow patterns without difficulties of principle. Introducing 
the stream function tp via U x = d y tp, U y = ~d x tp 7 the bulk 
equation in the liquid region, now nonlinear, reads [12]: 

% + T m = ^,T £ - + er^/fo- 1 ) 2 ^. (15) 

Here, the stream function tp is determined by tp^ + tp m = 
with the boundary conditions lim^oo (tp — tpo) = 0, 
where tpo = Pf(v — 1)£ is the Ivantsov-like solution for 
the stream function obtained with rf = 0, and 



1> e = -Pf{Zh)' 



(16) 



at the interface. The field equation in the solid region and 
the interface equations (3) remain unchanged. As before, 
we use (5) in the solid. Analogically to (4), we write 



(c\ + id v )tp = 0. 



(17) 



Equation (15), valid in the liquid region, does not factor- 
ize (not even asymptotically), but we may use Zauderer 's 
asymptotic decomposition method [13] to achieve a simi- 
lar reduction of order while keeping transcendentally small 
terms. In order not to overburden this first presentation, 
we linearize (15) about ipQ. Again, this is by no means a 
necessary step. We have performed the full analysis with- 
out this linearization, which leads to the same equations 
(22) and (23) given below. We first rewrite equations (15) 
and (17) together with the boundary conditions (3) and 
(16), and using (5), as a first-order system: 

w 6 + Awr, + e{Bw + Cw ? } = (18) 

with boundary condition 

aw + /3w £ = b , (19) 



e has been inserted into (18) as a bookkeeping variable, 
to keep track of the order of Zauderer's decomposition 
scheme. The formal procedure would be to write £ = 
—i + eX, r/ = 1 + eY and express the equations in terms of 
X and Y. In this presentation, we already decomposed the 
equations partially by writing Eq. (17). The parameter e 
is assumed small as <r < 1. It is possible to determine the 
scaling exponent for e in a more formal approach. A has 
the following eigenvectors: 



(20) 
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with eigenvalues i, i, 
tion of the form 



-i, respectively. We look for a solu- 



w = Mr n + ipr 12 + eNr 2 , 



(21) 



assuming e./V to be much smaller than M and ip. Ex- 
panding the resulting system of equations to first or- 
der in e, we obtain a set of decoupled first order equa- 
tions. Setting e = 1 and introducing characteristic co- 
ordinates s = —i(rj — 1), r = £ + i(rj — 1), we find 
tp = tp(s,r) = — Pf[Th(r)], and using the boundary con- 
ditions at the interface, we arrive at the two equations 

M = M(s, t) = l e T^ s2+rs ) J [(1 + tT)ft(r)]' 



- 2 [rh{r)\ 



i 1 — e 2 



(22) 



N 



S 

= N(s, f) = -^-f J ds' M(-s', 2s' + f) 



i^- I ds' [h{2s' + ?) + (2s' + f)h'(2s' + f )] 



+ l{o-K'(f)-[(l-if)h(f)}'} 



(23) 



where s = i(r) — 1), f = £ — i(rj — 1). Next, we require N 
to vanish for s — > ioo while keeping f fixed, f may be 
interpreted as the continuation of the variable £ into the 
lower half of the complex plane. Thus we write £ instead 
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of t, and after some manipulations Eq. (23) yields 



Studying the vicinity of £ = —i, we find 



ioo 



(i - it'K'hg) 

= an. (24) 



For Pf = 0, Eq. (24) reduces to (1). The calculation of 
the transcendental mismatch in the presence of a potential 
flow parallels the procedure for the flowlcss case. We first 
calculate the WKB solution of (24): 

h = B ie £(i - *o _f (i + cx p {""'5(0 + 

(25) 

with 5(0 given in Eq. (9). To obtain the inner equation, 
we follow Ben Amar [14] in her treatment of the flowless 
finite growth Pcclct number case and use the same scaling 
as in the case without convection, i.e., Eq. (10) with a = |. 
Equation (24) simplifies to 



tcf> + P 1 



J dt't'<j>{t') 



K . 



(26) 



with K defined in Eq. (12) and Pi = er?^-. To leading 
asymptotic order, 



P-t J dt't'<f>(t') ~ 2-1 1~ §, (i»l). 



(27) 



Linearizing equation (26) about (27) and performing a lo- 
cal asymptotic analysis near t — oo we obtain 



leading to 



= 2-H-i +A 1 (P 1 )g(t) (t-oo), (29) 



where Ai(Pi) is a nonlinear eigenvalue and a function of 
Pi. Matching (25) and (29), we obtain the tip slope 



2 V a 28 Im {^i} exp 1 - 



{^5(0)} 



(30) 

Obviously, the mismatch will remain nonzero except pos- 
sibly for isolated values of Pi , a case that may be excluded 
by numerical evaluation of Ai(Pi). Hence, there is no ad- 
missible solution. 

To include surface tension anisotropy, we have to replace 
k in (24) with Rk, where 



R = l-/?cos(40) = 1-/3 cos 



4 arctan 



Vs + tVs 



■ (31) 



t<jy + P 1 J At' t'cf>{t') =K + bHK . 



(32) 



where P = -2(0- l) 2 / (0 + l) 2 (</> + tf andfe = /3cr--. 

Demanding Im(Ai) to vanish, equation (32) constitutes 
an eigenvalue problem for the parameter b, to be solved 
numerically. Denoting the lowest eigenvalue by b = bo (Pi ) 
the selection criterion reads 



/(Pi)/? 7 / 4 , /(Pi) - MP 1 )]- J/4 



(33) 



We will not elaborate on the details of this solution, as 
the purpose of this article is only a demonstration of the 
method. Further results for the particular physical system 
will be discussed in a forthcoming paper. 

To summarize, we have introduced a method that com- 
bines matched asymptotics in the complex plane with the 
asymptotic decomposition of partial differential equations. 
This allows one to compute exponentially small terms be- 
yond all orders (of asymptotic expansions) for partial dif- 
ferential equations on free boundaries, as we have shown 
for dendritic growth in a forced potential flow. In com- 
parison with the BP approach [8] , ours has several advan- 
tages, the most important being that it paves the way for 
a rigorous nonlinear asymptotic analysis, which in some 
cases [9] seems to be the only one that gives even quali- 
tatively correct answers. We are not aware of any other 
method allowing this type of analysis with nonlinear field 
equations. Many problems to which the Kruskal-Segur 
method has been applied so far, including viscous finger- 
ing [15], the V 4 breather of particle physics [16], or capil- 
lary water waves [17], are actually free boundary problems. 
But for the method to be applicable, they first had to be 
recast as an ordinary differential equation or at least a sin- 
gle differential-integral equation. This is the reason, why 
only the simplest, in some cases even the most unrealis- 
tic, physical situations, have been studied by this method. 
Without the restriction, we anticipate a much wider class 
of problems to be tractable. Hence we expect our approach 
to open a new line of research into a plethora of hitherto 
untractable selection or solvability problems not only in 
crystal growth and similar problems of pattern forming in- 
terface dynamics (such as viscous fingering) , but in a vast 
number of situations reducible to free boundary problems. 
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